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ABSTRACT 

The  expanded  ensemble  Monte  Carlo  (EEMC)  simulation  method  has  been  applied  to 
calculation  of  the  chemical  potential  of  nanocolloidal  particles  in  the  presence  of  polymeric 
surface  modifiers.  Two  general  classes  of  surface  modifiers  have  been  studied  - nonadsorbing 
and  freely-adsorbing.  For  both  systems,  the  infinite  dilution  colloid  chemical  potential  was 
calculated  as  a function  of  the  colloid  diameter  and  the  modifier  chain  length.  The  colloid 
chemical  potential  was  found  to  decrease  with  increasing  modifier  chain  length  for  both  types  of 
modifiers,  albeit  for  different  reasons.  Empirical  power-law  scaling  relationships  were  found  to 
represent  the  simulation  results  well.  A physical  interpretation  was  proposed  for  the  power  law 
exponents  obtained  in  the  case  of  adsorbing  modifiers. 

INTRODUCTION 

Nanoscale  colloidal  dispersions  play  an  important  role  in  most  synthesis  and  processing 
steps  for  nanoscale  structural  materials.  Applications  include  catalysts,  adsorbents,  polymer 
nanocomposites  for  high-strength  barrier  and  flame  resistant  materials,  and  nanocrystals  for 
photonics  and  electronics.  These  nanostructures  can  be  prepared  not  only  via  nucleation  and 
arrested  growth  processes,  but  also  as  equilibrium  products,  allowing  thermodynamic  control  of 
crystal  lattice  structure  and  feature  size.  Many  experimental  approaches  utilize  short-chain 
organic  modifiers  to  promote  nanoparticle  stability  and  to  influence  the  lattice  structure  and  size. 
Examples  of  common  organically-modified  nanocolloids  include  noble  metal  nanoparticles 
stabilized  with  poly  (2-hydroxyl  methacrylates)  [1],  two-  and  three-dimensional  arrays  of  Pd  or 
CdSe  nanocrystals  modified  with  alkanethiols  [2,  3],  and  polymer-silica  nanocomposites  that 
mimic  natural  organic-inorganic  materials.  [4, 5]  Unfortunately,  the  present  knowledge  of 
organically-modified  nanoparticle  self-assembly  is  inadequate  for  precise  thermodynamic  control 
of  these  structures.  A breakthrough  in  rational  predictive  strategies  for  nanoscale  structural 
material  synthesis  would  be  the  development  of  molecular-models  that  account  for  the  effect  of 
surfactants  on  the  thermodynamic  and  transport  properties  of  nanoparticles. 

Towards  this  goal,  we  have  recently  developred  a novel  application  of  the  expanded 
ensemble  Monte  Carlo  (EEMC)  simulation  method,  which  allows  accurate  calculation  of  the 
chemical  potentials  of  organically-modified  nanoparticles.  Chemical  potential  (ji)  is  the  free 
energy  of  inserting  a molecule  (collection  of  molecules  in  the  case  of  a nanoparticle)  into  a 
system  of  interest  relative  to  a reference  state.  Knowledge  of  the  modifier  chain  length, 
concentration  and  particle  size  dependence  of  //  can  be  used  to  predict  self-assembly  and 
dispersability  (phase  behavior)  in  these  systems.  The  traditional  Widom  [6]  method  for 
calculating  fi  using  computer  simulations  works  well  for  small  molecules  and  moderate  densities, 
but  severe  difficulty  arises  for  polymers  and  other  molecular  aggregates  such  as  colloids,  due  to 
the  high  probability  of  steric  overlaps.  As  a result  the  Widom  method  samples  almost  exclusively 
the  high-energy  overlap  contributions  to  the  chemical  potential.  EEMC  is  a powerful  approach 


193 


that  overcomes  these  sampling  limitations  in  free  energy  calculations.  [7-14]Additional  ensemble 
variables  are  introduced  to  define  a reversible  and  efficient  path  between  two  desired  states  of 
free  energy  measurement.  These  states  are  sampled  according  to  probabilities  determined  from 
the  partition  function.  We  have  adapted  this  method  to  polymer-nanocolloid  mixtures  by 
inserting  nanoparticles  into  the  system  incrementally  in  a series  of  smaller  diameter  increments. 
[15] 

In  this  paper  we  have  applied  the  EEMC  method  for  calculating  nanoparticle  chemical 
potentials  in  the  presence  of  two  broad  classes  of  polymeric  modifiers:  nonadsorbing  and  freely 
adsorbing.  We  use  the  hard  sphere  (HS)  potential  to  model  nonadsorbing  modifiers  and  the 
Lennard-Jones  (LJ)  6-12  potential  to  model  freely  adsorbing  modifiers.  In  each  of  these  ca.ses  we 
calculate  infinite  dilution  nanoparticle  chemical  potential  as  a function  of  both  modifier  chain 
length  (n)  and  particle  diameter  (a,  ) and  discuss  the  physical  interpretation  of  these 
dependencies.  Finally,  we  present  empirically-derived  scaling  relationships  which  describe 
nanoparticle  chemical  potential  as  a power  law  function  of  n and  <7,  . 

THEORY 

To  adapt  the  EEMC  approach  to  nanoparticle-polymcr  mixtures,  the  particle  diameter  (the 
“expansion  variable”)  is  allowed  to  vary  between  zero  (fully  removed)  and  full  size  via 
increments  and  decrements.  These  increments  or  decrements  are  accepted  according  to  a 
Metropolis  type  acceptance  criterion.  A predefined  number  of  intermediate  diameters  of  the 
particle  are  allowed  by  prescribing  them  as  the  different  “states”  of  the  EE.  Incremental 
chemical  potentials,  calculated  at  each  diameter  describe  the  free  energy  of  increasing  diameter 
from  di  to  The  full  chemical  potential  is  recovered  by  summing  the  incremental  values  from 
zero  to  the  desired  particle  size.  The  statistical  mechanical  details  of  this  method  have  been 
presented  in  a recent  publication.  [15] 

SIMULATION  DETAILS 

The  simulations  were  performed  in  a cubic  box  with  periodic  boundary  conditions.  In  all 
cases  our  system  consisted  of  a single  colloidal  particle  dispersed  in  a dilute  polymer  solution. 
Polymer  chains  are  fully  flexible  and  consist  of  tangent  segments  with  a constant  bond  length  of 
(7p  (bead  spring  model).  The  polymer  segment  density  Pp  is  maintained  constant  for  all  the 
simulations  at  Pp  = 0.05,  which  corresponds  to  a packing  fraction  of  <t>  = 0.026 . The  chain  length 
of  the  polymer  was  varied  from  n = 2 to  a maximum  of  « = 60.  The  effect  of  colloid  size  on  the 
chemical  potential  was  also  studied  by  varying  a,  from  /q,  to  I5ap  . In  all  of  the  simulations, 
the  box  length  L (typically  40oj,)  was  greater  than  ;i/2  and  ct  . to  prevent  artificial  intrachain  or 
intraparticle  interactions. 

RESULTS 

Nonadsorbing  modifiers  (HS  potential) 

The  HS  potential,  used  here  to  model  nanoparticles  in  the  presence  of  nonadsorbing 
modifiers  in  a good  solvent,  gives  purely  the  entropic  or  configurational  contribution  to  the 
chemical  potential.  The  HS  model  system  provides  a stringent  test  of  EEMC  chemical  potential 
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calculation  since  the  configurational  contribution  is  the  most  difficult  part  of  the  calculation.  In 
addition  it  has  been  shown  that  self-assembly  can  be  driven  by  entropic  depletion  forces  in  the 
presence  of  nonadsorbing  polymer.  [16, 17]  Rgure  1 shows  the  infinite  dilution  chemical 
potential,  versus  the  colloid  diameter  ac  for  polymer  chain  lengths  of  5 < n < 60.  In 
general  increases  as  a cubic  polynomial  in  for  all  chain  lengths,  reflecting  the  increase  in 
the  probability  of  overlaps  and  excluded  volume.  Figure  1 compares  the  simulation  results  with 
the  predictions  of  the  PRISM-based  colloid-polymer  model  of  Fuchs  and  Schweizer  (FS  model). 
[18,  19]  The  FS  model  compares  very  well  with  the  simulations  except  that  at  chain  lengths 
above  30  the  FS  model  slightly  overpredicts  the  chemical  potential,  and  at  chain  lengths  less  than 
30,  the  predictions  are  slightly  lower  than  the  simulated  results.  Minor  quantitative  deviations 
are  due  to  the  assumed  Gaussian  intramolecular  polymer  structure  factor  in  the  model.  Results 
from  a field-theoretic  (FT)  model  are  also  shown  for  comparison.[20] 

Figure  2 shows  the  effect  of  the  polymer  chain  length  on  the  infinite  dilution  chemical 
potential,  Pfj/' , for  different  colloid  diameters  Oc.  The  chemical  potential  decreases  with  chain 
length  for  all  particle  sizes.  For  a constant  polymer  segment  density,  longer  chain  length  (i.e., 
lesser  number  of  chains)  results  in  less  excluded  volume  for  the  particle  due  to  chain 
connectivity  con.straints  on  the  location  of  polymer  segments.  This  explains  the  decrease  of 
chemical  potential  with  chain  length.  From  Figure  2 we  observe  that  for  all  particle  sizes,  there  is 
a sharp  decrease  in  colloid  chemical  potential  from  « = 2 to  n = 20.  In  this  regime  inserting  the 
colloid  particle  requires  the  free  energy  of  creating  additional  volume  by  doing  work  against  the 
polymer  fluid  osmotic  pressure.  However  for  larger  chain  lengths  only  an  internal  section  of  the 
polymer  chain  needs  to  be  rearranged  as  the  colloid  particle  sees  only  a part  of  the  chain.  Hence 
we  see  that  for  n > 20  the  change  in  chemical  potential  is  comparatively  small.  This  is  an 
interesting  result  since  it  tells  us  that  firstly,  even  in  the  absence  of  attraction,  a large  reduction  in 


Figure  1.  Infinite  dilution  HS  colloid 
chemical  potential  vs.  colloid  diameter  in 
a dilute  HS  polymer  solution  at  different 
polymer  chain  lengths  (n).  Dotted  line: 
FS  model.  Solid  line:  FT  model. 
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Figure  2 Infinite  dilution  HS  colloid 
chemical  potential  vs.  polymer  chain 
length  in  a dilute  HS  polymer  solution  for 
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the  colloid  chemical  potential  can  be  brought  about  by  adding  short  chain  polymers  (modifiers) 
and  secondly,  increasing  polymer  chain  length  beyond  a particular  value  yields  little  additional 
reduction  in  chemical  potential.  In  addition,  the  ehemical  potential  reduction  effect  from  adding 
polymer  chains  increases  as  particle  size  increases. 

Freely  adsorbing  modifiers  (L.I  6-12  potential) 

The  U 6-12  potential,  used  to  model  freely  adsorbing  homopolymer  modifiers  incorporates 
attraction  and  chain  adsorption  effects.  The  LJ  model  captures  the  essential  physical  features  of 
many  organic  molecules  and  remains  sufficiently  general  to  be  of  broad,  predictive  value.  All  the 
simulations  were  performed  at  a reduced  temperature  of  T*  = 3.0  (T*  = Iks/s)  and  U energy 
parameters  Eir  = Epp  = £cp  = 1,  where  ^c,£,,,,and  f^,,are  the  colloid,  polymer  and  colloid-polymer 
LJ  interaction  parameters  respectively.  Figure  3 shows  a plot  of  infinite  dilution  chemical 
potential,  Pfi/'  versus  the  colloid  diameter  o,  for  polymer  chain  lengths  of  5 < n<  30.  For  all 
chain  lengths,  decreases  with  o-  and  shows  a cubic  dependence  on  Oc-  The  chemical 
potential  decreases  with  o;  because  the  LI  attractive  energy  increases  linearly  with  o,  and  this 
results  in  stronger  attractive  (negative)  interactions  between  the  particle  and  the  polymeric 
modifiers  and  thereby,  a lower  chemical  potential.  These  results  demonstrate  that  in  the  case  of 
strongly  adsorbing  modifiers,  attractive  energetics  dominate  the  repulsive  and  excluded  volume 
interactions  that  are  predominant  in  nonadsorbing  modifier  systems. 

The  effect  of  chain  length  of  adsorbing  modifiers  on  particle  chemical  potential  is  shown  in 
Figure  4.  As  in  the  case  of  nonadsorbing  modifiers,  here  too  we  observe  that  for  all  particle 
diameters,  decreases  with  increasing  chain  length.  However,  more  than  the  decrease  in 
excluded  volume  due  to  chain  connectivity  constraints,  it  is  chain  adsorption  on  the  particle 
surface  that  leads  to  lowering  of  particle  chemical  potential  in  this  case.  As  we  increase  the  chain 
length,  the  entropy  loss  associated  with  chain  adsorption  decreases  and  hence  more  chain 
segments  get  adsorbed  onto  the  particle  surface.  This  results  in  the  particle  getting  a larger 
number  of  attractive  negative  interactions  and  thus,  a lower  chemical  potential. 

Scaling  relationships 

The  chemical-physics  of  nanocolloid-polymer  mixtures  can  be  described  by  simplified 
scaling  relationships  which  capture  the  dependence  of  thermodynamic  properties  on  physical 
molecular  parameters.  These  relationships  reflect  the  physics  that  is  observed  in  the  simulations 
and  could  be  used  to  compare,  interpolate  of  extrapolate  experimental  data.  With  this  motivation 
we  performed  power  law  regression  of  the  chemical  potential  data  for  the  two  classes  of 
modifiers  using  a,-  and  n as  the  functional  molecular  parameters  and  arrived  at  the  following 
relationships: 

Nonadsorbing  modifiers  Pjtc'’'=  0.13Oc^'^n'”‘”* 

Freely  adsorbing  modifiers  PMc'*  = - 0.012Oc^  '”n"  '''' 

For  nonadsorbing  modifiers  modeled  with  the  HS  potential  it  has  been  shown  that  the 
chemical  potential  shows  a cubic  polynomial  dependence  on  cr,  !Rg  where  Rg  is  the  radius  of 
gyration  of  the  polymer.  [15,  18,  19]  Hence  the  empirical  power  law  scaling  exponents  shown 
above  have  little  physical  significance  and  are  a mathematical  con.sequence  of  reducing  a cubic 
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polynomial  to  a power  law  function  of  the  independent  molecular  parameters.  However,  for  the 
freely  adsorbing  modifiers,  we  believe  that  a physical  import  can  be  attributed  to  the  exponents 
in  the  scaling  relationship.  The  cubic  dependence  on  CTc  is  essentially  a product  of  a linear  term 
(a;)  and  a square  term  (o;  ^).  The  linear  dependence  represents  unit  particle-polymer  interaction 
energy  since  the  LJ  attractive  energy  varies  linearly  with  the  particle  diameter.  The  square  term 
is  a surface  area  effect  which  is  indicative  of  the  number  of  chains  that  adsorb  on  the  particle  i.e., 
the  number  of  particle-polymer  interactions.  The  chain  length  exponent  arises  out  chain 
adsorption  and  represents  the  thickness  of  the  adsorbed  layer.  For  chains  that  are  anchored  on  a 
surface  in  a theta  solvent,  the  mean  field  theory  predicts  that  the  thickness  of  the  adsorbed  layer 
scales  with  [21]  In  the  case  of  freely  adsorbing  modifiers,  the  chains  would  tend  to  swell  a 
little  more  compared  to  tethered  chains  and  hence  one  would  expect  that  the  exponent  would  be 
slightly  lower  than  0.5.  hi  addition,  since  the  adsorbed  chains  are  closely  packed  near  the  particle 
surface,  the  attractive  and  repulsive  interactions  between  the  chain  segments  would  be  relatively 
balanced,  thereby  simulating  close  to  theta  solvent  conditions. 

CONCLUSIONS 

The  EEMC  method  was  adapted  and  applied  to  calculate  the  infinite  dilution  nanopaiticle 
chemical  potential  in  the  presence  of  organic  surface  modifiers  as  a function  of  particle  diameter 
and  modifier  chain  length.  Two  broad  classes  of  organic  modifiers  were  explored:  nonadsorbing 
modifiers  and  freely-adsorbing  modifiers.  For  both  cases  it  was  shown  that  the  addition  of  short 
chain  organic  modifiers  (n  < 20)  brings  about  a substantial  reduction  in  the  chemical  potential, 
compared  to  pure  monomer  solvent  at  the  same  volume  fraction.  The  dependence  of  particle 
chemical  potential  on  polymer  chain  length  and  particle  diameter  wm  found  to  be  represented 
well  by  empirical  power  law  scaling  relationships  of  the  form  a afit  for  both  types  of 
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modifiers.  For  nonadsorbing  modifiers  the  Oi  and  n dependencies  arise  purely  out  of  excluded 
volume  and  entropic  contributions.  However  for  adsorbing  modifiers,  the  LJ  attraction  energy, 
surface  area  and  chain  adsorption  are  responsible  for  the  dependence  of  ' on  (T,  and  n.  These 
scaling  relationships,  though  empirical,  could  be  used  to  guide  experiments  and  efforts  towards 
theoretical  development  of  chemically  specific  equation.s  of  state  for  these  complex  nanoparticlc- 
polymer  mixtures. 
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